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ABSTRACT 

The network design problems have been widely investi- 
gated m the mathematical programming literature in the last 
decade. In this thesis, we study one such design problem known 
as Optimum Communication Spanning Tree (OCST) problem. We 
develop an exact algorithm based on branch and bound approach 
and heuristic algorithms for solving the OCST problem. We first 
derive two new lower planes for the network design problem 
through combinatorial arguments. One of the lower planes is 
used in computing lower bounds of subproblems in the branch and 
bound algorithm. Reoptimization of subproblems is extensively 
used to reduce computations. The heuristic algorithm is a two- 
phase algorithms tree-building algorithm and tree-inprovement 
algorithm. The complexity of these algorithms is analysed. We 
also investigate a heuristic algorithm which is a branch and 
bound algorithm with inexact lower bounds. 

All these algorithms are programmed in FORTRAN IV and 
implemented on DEC-10 computer system. The algorithms are 
tested on randomly generated Euclidean and non-Euclidean 
network problems. The branch and bound algorithm is able to 
solve 40 node, 69 arc Euclidean network problem m about 37 
minutes of execution time. Non-Euclidean problems are solved 
more easily. The two-phase heuristic algorithm has produced 
excellent results. All the problems tested are solved 



optimally. T he computational times are also very reasonable . 
Problems of size 100 nodes and 1,000 arcs are solved in about 
30 seconds of execution time . The performance of the heuristic 
algorithm based on inexact bounds is not encouraging. 



CHAPTER I 


INTRODUCTION 


1.1 INTRODUCTION: 

In recent years, much interest has been focussed on 
various problems that arise in the area of network design. In 
its simplest form, the network design problem may be stated as 
the building of a connected sub-network, by selection from 
the set of arcs of a network, that minimizes the weighted 
sum of shortest path distances between all pairs of nodes of 
the sub-network, subjected to a budget constraint which limits tlr 
number of arcs that may be included in the sub-network. Such 
problems find application in transportation planning, communi- 
cation system planning, water resource planning and distribution 
planning. 

In this dissertation, we consider one such network 
design problem known as Optimum Communication Spanning Tree 
(OCST) problem. Given an undirected network and the amount of 
communication between each pair of nodes, the OCST problem is 
to determine a spanning tree for which the total cost of 
communication between all pairs of nodes is minimum among all 
spanning trees. Such problems are encountered very frequently 
in designing road, communication and computer networks. In 
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coimnunication networks, the nodes represent cities and the 
amount of communication may be taken as the number of telephone 
calls between pairs of cities. The total cost of communication 
decreases if more number of arcs are included in the sub- network 
but the construction cost limits the inclusion of more number 
of arcs. Thus a sub-network is to be built which barely connects 
all the cities and at the same time the total communication cost 
is minimum. This results in building a spanning tree that 
minimizes the total cost of communication among all spanning 
trees. 

In all these applications, the actual design problem is 
usually subjected to many constraints but the solution of the 
OCST problem may be used as a measuring standard for the effi- 
ciency of proposed designs. The OCST problem has been shown to 
be HP-complete [l4] even for the special case when the amount 
of communication is same for all pairs of nodes. 

In this dissertation, we develop exact and heuristic 
algorithms for the OCST problem. The exact algorithm is based 
on the branch and bound approach for solving combinatorial 
optimization problems. A new lower bounding strategy is 
developed for the general network design problem which is used 
in the branch and bound algorithm. The heuristic algorithm is 
a two-phase algorithm - the tree-building algorithm and the 
tree-improvement algorithm. The tree-building algorithm starts 
with a seed node and successively builds the tree by adding one 
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arc each time such that a good feasible solution is obtained. 
The spanning tree thus constructed is fed to the tree improve- 
ment algorithm which examines each tree arc one by one and 
tries to interchange it with an appropriate non-tree arc so 
that the cost of communication decreases. The algorithm ter- 
minates when no possible interchange yields improvement. 
Another heuristic algorithm, named Global heuristic, is also 
developed. It is essentially the branch and bound algorithm 
with the modification that in-exact lower bound is used. This 
is an adaptation of the heuristic proposed by Dionne and 
Florian [5] for the network design problem and was claimed to 
be giving excellent results. 

All the suggested algorithms have been programmed and 
tested extensively on randomly generated network problems. The 
branch and bound algorithm is able to solve 40 nodes, 69 arcs 
problem. For all the problems solved by the branch and bound 
algorithm, the two-phase heuristic algorithm gave the optimal 
solution. The global heuristic algorithm was observed to be 
taking more time when compared to the two-phase heuristic 
algorithm. Detailed computational results are presented. 

1.2 PRELII iINARIESs 

Some notations and well-known concepts of graph theory, 
used throughout the thesis, are given below. 

A g raph G = (N, A) consists of a finite set N of 
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elements, called nodes, and a set A of pairs of nodes called 
arcs,. An undirected gra ph is a graph in which the arcs have 
no direction. 

A graph in which numbers are associated with nodes and 
arcs is called a network. A path in G a (N, A) is a sequence 
il,i 2 ,..., i r of distinct nodes of N such that (i^, i^ + ^ ) e A> 
for each k = 1 ,..., r- 1 . A cycle, is a path together with an 
arc (i r , ij) . 

A graph G’ = (N 1 , A 1 ) is a subgraph of G = (N, A) if 
N* <p, N and A 1 c: A. A graph is connected if there is a path between 
any two nodes of the graph. A tree, of k nodes is a graph which 
is connected and has k - 1 arcs. A sub-graph which is a tree 
and contains all vertices of a graph is called a ^anning^tree,. 

(In a spanning tree, there is a unique path between any two nodes) 
Let, 


m = 

|A|, n = |N| 


r i 3 

amount of communication between the pair of nodes i 


and 


d « ■ 

length of the arc [i,o], i £ 3 ( = “ if there 

no arc between 1 and 3 ) 

is 

d n - 

0, for i = 1, . . . , n. 


x = 

a subset of A 



= shortest distance between the pair of nodes 

and 3 over tin. graph G(N, %) , 

i 
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The following algorithms are extensively used in the 
present work and are given below. 

1 , Floyd's [ 7 ] Algorithm for Finding Shortest Path Distances 
between all pairs of nodes: 

Procedure FLOYD (n, d, u); 


begin 


Set u := d . 


il H 7 

for h;= 1 to n do 


for 1 to n-1 do 
for 3*= i+l to n do 

if u i0 > U ih + u hj then 
begin 


U i3 : “ U ih + U ho ? u : 


¥ i: = 1,..., n, and ¥ 1,..., n. 


2 . Dionne and Florian’s [ 5 j Algorithm for Updating the 
Shortest path distances when an arc length changes: 

Let [p,q] be the modified arc, its length being d 

P 

before modification and d^ afterwards. The algorithm 
obtains the modified u matrix when the length of the arc 
changes from to d^ procedure ADDARC for the case when 
dpq > and procedure DELETEARC for the case when d < d^. 
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Procedure ADDARC j 
begin 

Nl; = 0 ; N 2 s= 0 ; 
for 1 to n do 

lf U iq > u ip + d pq. then 
begin 

Nl;= Nl u UK u iq := u lp H. <^ q ; u ql ••= u lq 
end; 

for 32= 1 to n do 

t 

if u , > d „ + then 

pa pq qa 

begin 

H2: = N2 U UK d' q + u q;)i u jp; « u p;) 

end; 

for every (1,3) e Nl x N 2 do 

if u. _ > u_ + u then 
10 ip po 

begin 

u s = u,^ + ti ; u. 2= u. 

13 ip pa 0^* ia 

end 


end; 

Procedure DELETEARC; 
begin 

Nl? = 0 ; N 2 S= 0 j Rs= 0 ; 
for i.= 1 to n do 

if u iq = u ip + d pq then Ml:= Nl U {i}; 
for o?= 1 to n do 


if Upj = d pq + u q;J then N2:= N 2 U{ 3 }; 
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for every (i,j) e N1 x N2 do 

if - u i P + Si + % then 

begin 


R:= r u {(i,d)b Ujlj*" d 13 i u i 3 ; 

end; 


V s= d' id i« d ;u s = d iu ~ 
pq» qp pq» pq pq» qp 

for hs= 1 to n do 

for every (i,j) e R do 

if u ih + u hj < u ia then 
begin 

u ij ? = U ih + u ho»* u 3 i S “ u i3 
end 



end; 

This procedure DELETEARC is modified to delete a number 
of arcs at a time and its given as procedure DELETEARCS . 


Procedure DELETEARCS (P,q)| 

{Modifies the shortest distance matrix u when all arcs 
[p,q], P e P are removed from the network, R is the 
set of retained pairs} 
begin 

R— 0? 

for every p e P do 
begin 

NlS as 0j N2!= 0J 

for i!= 1 to n do if u, * u. + d then 

iq. rp pq 

Nls= N1 U {i}; 
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for o!= 1 to n do if then 

N2s « N2 U {j}; 

for every (i,o) e HI x N2 do 
begin 

if u ip + + s u i 0 then R * = R u 

end; 

Sq 5 " d pq*‘ d qp 5= d pq 

end; 

for every (i,d) e R do 
begin 

u ij : = V u ji != u i3 

end; 

for hs= 1 to n do 

for every (i,i) e R do 

if u ih + "ha < u ia then 
begin 


u ij ; = u ih + u h D '* u ox i= U i0 


end; 


end; 


3. Kruskal’s [15] Algorithm for Minimum Spanning Tree Problems 


Procedure MST; 
begin 

T?= 0; ES = A; 

While |Tl < n-1 and E ^ 0 do 
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begin 

find an edge e in E with minimum weight) 

E: = E— { e } j 

if T U {e} has no cycle then Tu T U {e} 
end 

end) 

The maximum spanning tree problem is solved using 
Kruskal's algorithm and the data structures given in [24] 
are employed for implementing this algorithm. 


1.3 STATEMENT OF THE PROBLEM: 


cimum C ommuni c at ion 


>annu 


Given an undirected network and the amount of communica- 
tion, hy , between each pair of nodes i and j, the OCST problem 
is to determine a spanning tree, for which the total cost of 
communication between all pairs of nodes is minimum among all 
spanning trees. The cost of a given spanning tree is defined 
as follows. For a pair of nodes i and 3 , there is a unique 
path in the spanning tree between 1 and The length of the 
path is the sum of lengths of arcs in the path. The cost of 
communication for the pair of nodes i and j is r^ multiplied 
by the length of the path. Summing over all n(n~l )/2 pairs 
of nodes, we get the cost of the spanning tree. This can be 
stated as, 
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-i XX XX 

Min. 4 I I r. ..u. . (A) 

* i=l 3=1 3 3 

A is a spanning tree. 

Two special cases of the OCST problem may he stated ass 

1. Optimum Requirement Spanning Trees The arc lengths are 
all equal to 1, while the requirements r^ are arbitrary. 

2 . Optimum Distance Spanning Trees The r . - are all equal to 1, 

X J 

while the are arbitrary. 

The general OCST problem has been formulated by Hu [l3]. 
He has also developed an O(n^) algorithm for solving the opti- 
mum requirement spanning tree problem. The optimum distance 
spanning tree problem has been shown to be HP- complete by 
Johnson, Lenstra and Rinnooy Kan [l4] by showing that the 
Exact 3-Cover is reducible to the above problem, 

1.4 OUTLINE OF THE THESIS s 

A brief chapter by chapter outline of the thesis is given 
in this section. 

In Chapter XI, existing lower planes for the network 
design problem are described and two new lower planes are 
proposed. An 0(mn ) lower plane by Hoang [ll] and an 0(n H ) 
lower plane by Gallo [9] have been proposed for the network 
design problem in the literature. A new lower plane of worst 
case complexity O(n^) which is tighter than that of Gallo is 
developed in Section 2.3- Another lower plane of O(n^) 
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complexity is proposed in Section 2.4. A comparison of lower 
bounds by all the methods for randomly generated problems is 
presented in Section 2.5. 

In Chapter III, we develop an exact algorithm based on 
the branch and bound approach and two heuristic algorithms - 

(i) Two-Phase heuristic and (ii) Global heuristic* The new 
lower plane developed in Section 2.3* is used for lower bound 
calculation. Some of the important features of the branch and 
bound algorithm are: 

(i) Sensitivity analysis is used to update the shortest path 
distance matrix whenever an arc length changes (increases or 
decreases). 

(ii) Reoptimization is done to reduce the computations subs- 
tantially. 

(iii) The new lower plane of Section 2.3 along with the maximum 
spanning tree algorithm is used in the lower bound calculation. 

(iv) The arc selected for branching is such that the pruning 
is more effective. 

(v) Dejth first search strategy is employed to reduce the 
storage requirements. 

Detailed discussion of all these aspects is presented. A nume- 
rical example is given to clarify the above steps. 

In the first phase of the Two-Phase heuristic algorithm, 
a spanning tree is constructed by adding one arc each time such 
that a good feasible solution is obtained. In the second phase. 
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each tree-arc is examined one-by-one and an interchange of the 
tree-arc with an appropriate nonr-tree arc is made if the cost 
of communication decreases by such an interchange. The algorithm 
terminates when no such interchange yields improvement . The 
global heuristic algorithm is same as the branch and bound algo- 
rithm except that an in-exact lower bounding is employed. Because 
of the increased lower bound, the algorithm takes less computa- 
tional time as compared with the branch and bound algorithm, 
but the solution may not be optimal. 

In the same chapter, we present the computational perfor- 
mance of branch and bound, two-phase heuristic and global 
heuristic algorithms. The results are found to be quite encou- 
raging. The branch and bound algorithm is able to solve 40 nodes 
and 69 arcs problem in 37 min. of CPU time. For all the problems 
solved by the branch and bound algorithm, the two-phase heuris- 
tic algorithm gave the optimal solution. The global heuristic 
algorithm was observed to be taking excessive times when com- 
pared with the two-phase heuristic algorithm. Computer programs 
for all these algorithms are written in FORTRAN - IV and 
implemented on DEC-1 090 computer system. Detailed computational 
results are presented. 



CHAPTER II 


NEW LOWER PLANES FOR THE NETWORK DESIGN PROBLEM 

2.1 INTRODUCTIONS 

The network design problems have been widely studied in 
literature [20] . These are difficult combinatorial optimiza- 
tion problems and have been shown to be NP-hard [14], Benders 
decomposition [4,6,10,12,18,19,22] and branch and bound appro- 
aches [1,2,3,5,8,11,16,17,23] have been the two most effective 
optimization based solution techniques. However, despite these 
efforts only moderately sized problems (upto 30 nodes and 60 arcs 
have been solved within reasonable computational times. In 
this chapter, we develop two new lower bounding methods for 
the network design problem, whose use in the branch and bound 
algorithms might enable us to solve larger sized problems. 

The network design problem, addressed in this chapter, 
is to identify a subnetwork of a given network G satisfying a 
budget constraint, over which the cost of transportation is 
minimum among all such subnetworks. A subproblem in the branch 
and bound algorithm for the network design problem consists of 
a set of arcs, I, to be included in the network, another set 
of arcs, E, to be excluded from the network, and rest of the 
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arcs, U, unassigned. A naive lower bound jj, can be obtained 
by determining all pair shortest path distances u. oven, the 
network consisting of the arcs in I U U and letting 




o 


1 

2 


n n 

i=l o-l 3 J 


where r^^ is the amount of interaction between the pair of nodes 
i and ;j. However, the budget will not permit all unassigned 
arcs to be included in the network* Therefore, some arcs in U 
will have to be deleted leading to increase in the shortest path 
distances and, accordingly, the lower bound value. Hoang [ll] 
has suggested an 0(mn ) method which considers the effect of 
possible arc deletion on m pairs of nodes* Recently, Gallo [9] 
has proposed an O(n^) method which considers the effect of arc 
deletion on n(n-l)/2 pairs of nodes. The first lower bounding 
method suggested by us can consider the effect of arc deletion 
on each pair of nodes more than once and, consequently, obtains 
a bound which is always sharper than that of Gallo. The com- 
plexity of this method^is shown to be O(h^). The second lower 
bounding method suggested by us is of complexity O(n^) and produ- 
ces a bound which is generally worse than that of Gallo, but 
much sharper than that of Hoang. Computational results are 
also presented* 


2.2 NOTATIONS? 


We use the notations similar to those of Gallo [9] . For 
the sake of completeness, we give them below. Let, 
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G = (N, A) be an undirected network, 

n = |N| , 

m = |A|, 

X = a subset of arc set A, 

= the k-th arc, also represented as 
c k “ cos "^ constructing the arc a^, 
b = budget for constructing arcs, 

« length of the arc a^., 

r. = amount of interaction between the pair of 

nodes i and j, 

Q = {(i,a)s 1 < i < 0 < n}, the set of all pairs of 

nodes , 

u ij(^)“ distance between the nodes l and j over 

the graph C? = (N, X), 

g. (X)= the second last node on the shortest path in 

J.J 

G(N, X) from node l to node q if i / 3, and 
0 if i = j. These indices can be obtained while 
computing all pair shortest path distances, 

X = a subset of X. Define a characteristic vector of 
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The following results can be easily deduced from the 
above definitions! 


Pronosition i s If X c X , then n(X , (1,3)) < rc(X , (i,j)), 

¥ (i, 3 ) e a. 

Pronosition 2 s L (X - X, Q) *= p + 2 r. . .n(X, (i, J) ) . 

(i,0)eQ XD 

The Network Design Problem (NDP) may be defined as to 
select a set X c A such that 


2,3 


(i) L(A - X, Q) is minimum, and 

m 

(ii) x e S = {xi cx > 2 c, - b, and A - X is connected}, 

~ k=l K 


LONER PLANES 


A subproblem in the branch and bound, algorithm for the 
NDP consists of solving the problem over G = (N, X), where 3C c A. 

A linear function h(x) = p Q + p-x defined over if 1 is a lower plane 
for NDP if h(x) <_ L(%-X, Q) for every x £ S. Given any two lower 
planes h(x) and g(x), we say that h(x) dominates g(x) if h(x) g(x) 
for every x e S. 


2.3*1 Existing Lower Planes 

Hoang’s [ll] lower planes h^Cx) = p Q + p^x is a lower plane 
for the NDP, where 


“k - r i k 3 k - 1,(a k’ (i k-Jk )) > V 

o 

Gallo’s [93 lower planes h (x) « p Q 
for the NDP, where 


1c X 9 0 » 9 f nu 

2 

+ \x x is a lower plane 


( i 1 0 / £ 
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and, 

Q k = {(i,3) e Q* Tt(a k ,(i,j)) > 0 and te^, (i,3)) “ °» 

¥ p = 1»..«, k-l}, ¥ k *» m. 

P 1 T 

Clearly, h (x) dominates h (x) because if > 0 for 

some k, then Q k always includes the node pair (i k , 3 k ). Intuitively, 

h^(x) considers the effect of arc deletion on m pairs of nodes, 

whereas h^(x) considers n(n-l)/2 such interactions. Accordingly, 

Gallo 1 s numerical investigations have indicated that h^Cx) is 

much sharper than h^Cx) , 


2.3.2 The Lower Plane LPls 

We derive another lower plane for the NDP which considers 
the effect of arc deletion on each pair of nodes more than once 
but, in effect, it considers n(n-l)/2 interactions with highest 
magnitude. For each (i,a) e Q, define 


Hi ■ °’ 


and. 


3 k 
3 10 


_k 


*= max. {u(a , (i,d))}, ¥ k = 1,..., m, 
l<p<k p 

a max. {0, % (a. , (i,o)) - ¥ k = 1,..., m. 


mew*. # >w 


The following result easily follows from the above definitions. 

k k 

\% For any arc a^, if w^ > 0, then Z w^^ « 7i(a k , (i,j 

Theorem i s h^(x) = |i Q + ^ x is a lower plane for the NDP, where 

pj = Z r,_ wf ., ¥ k = 1,..., m. 

K (i,3)ea 13 10 
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Proof; Consider any X such that its characteristic vector x e S. 
Further, consider any node pair (i,o)» Let s be the highest 
index such that x_ = 1 and w^. > 0. If no index satisfies 
these conditions, then the weight contributed by this node pair 


to deleted arcs is zero. Clearly 


m k 
1 w ij * x k 


s 

E 


k=l 


k k 

w ii* x k — ^ w ii* 

k=l 10 K k=l 10 


Using Proposition 3 and Proposition 1 we get 

m lr s k 

2 x, < E wf = 7t(a_, (i,d) < Tt(X, (i,a)). 

k=l 1J K k=l 13 s 

Multiplying this by r^ and adding for all (i,j) e Q, 
we obtain, 


E 

(i,d)eQ 


,r ia 


m 

( I 

k=l 


"ij 




E r. . ,tc(X, (i, o) ), 

(i,j)eQ 


which after rearranging terms and using Proposition 2 becomes 

m v „ 

I ( Z r ±1 wj ) x k < L (A - X, Q) - VI , 
k=l (i,d)sQ 10 13 K ° 

or 

x <_ L (A - X, Q), 

Theorem 2t The lower plane h^(x) dominates the lower plane h^(x). 

Proof; Every element (i,a) e has uCa^, (i,j)) > 0 and 
7t(a , (i,a)) =0, V p = 1,..., k-1. Hence = 0 and 

•n. ** ^2 

« it(a k , (i,o)). It immediately follows that p. ^ 

V k = 1,..., m. Further, there may be some (i,j) e Cl such that 
> 0 and 7c(a k , (i,j)) > which will make the inequality 

^ k k strict * 
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A procedure, named LP1, to determine h^(x) is stated 

"below. The procedure deletes each arc a^ e % one by one, retains 

the pairs of nodes which might be affected by deleting the arcs 

in a set S, applies Floyd's algorithm [7] for each retained 

pair to determine the length of the changed shortest path between 

3 

the pair, and accordingly determines the weight to be asso- 
ciated with the arc a^. 

Pro-position 4; The computational complexity of determining 
h 3 (x) is 0(n 4 ). 

Proof: The shortest path between any pair of nodes consists of 
atmost n - 1 arcs and deletion of only these arcs can change the 
shortest path between the pair. Hence each pair of nodes will 
be retained atmost n - 1 times while deleting arcs one by one. 
Since there are n(n-l)/2 pairs of nodes, the total number of 

Q 

retained pairs will be atmost n(n-l) /2. 

The all pair shortest path lengths can be obtained in 
O(n^) operations. The operations in loop-1 are of 0(n 2 ) and 

p 

since it is executed atmost m times, its complexity is O(mn^). 

Ihe complexity of loop-2 is O(n^), as 0(n) computations are per- 

formed for every retained pair which are atmost n(n-l) /2 in 

number. The complexity of loop-3 is clearly O(n^). Hence the 

/ 4n 

overall complexity of the procedure LP1 is 0(n ). 



procedure LPlj 
■begin 

compute all pair shortest path lengths 1 s and the 
corresponding ' s; 

Set p..t = 0, V i» = 1,,.., n, ¥ j:= 1,..., n; 

n n 

compute 1 S r...u ; 

° „ i=l o=l 10 

for each a k e A do 

begin 

p: = i k j qs = j k ; 

if V ■ d pq then 
begin 

temp: = d^, d^s = dqp , . d^; 

Nls = 0,* N2s = 0j Ss = 0; 

for is = 1 to n do if = p then Nls = N1 U {i}; 
for o ! = 1 to n do if g_. p = q then N2t = N2 U (o}» 
for every (i,o) e N1 x N2 do 
if u. . = u. + temp + u then 

il ip qo 

begin 


S: = S U { ( i, o ) } i 

V. -S = u. I V . - s = V. . j 
10 . 10 ’ ,01 10 1 
end; (loop - 1} 

for hs = 1 to n do 

for every (i,o) e S do 

lf u ±h + u h3 < u ld then 


u. 


10 


= d. 


id’ 0i 


u . s = u. 


id 


begin 


= u ih + u hd 5 


U i0 

end; (loop - 2j 


u n ! 


= u. 


10 


weights = 0; 

for every (i,d) £ s do 

^ > p±d then 

begin 

weights = weight + r iD * ^ u io“ v id”^id ^ * 


* u.._-v,^; p^S 


end 

end 


enilO 10 10 

u id; = v id; u o 

end; {loop-3} 




= u 


i0 


P k * = weight; d pq s = temp; d^ ;= d. 


p<i 


end; 
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2* 3. 3 The Lower Plane LP2 s 

¥e now derive another lower plane for the NDP vdth 
O(n^) worst case complexity. The following additional numbers 
are defined for each (i,j) e CU 

P. . = the set of arcs in the shortest path between 

the nodes i and (i f o), 

p i,j = ^ P io^ 

v^ = the second shortest path length between the 
nodes i and 

A ij = v io “ u id’ 
v f 1. if P. .. contains the arc a, 

ah 4 13 k (1) 

0 t °» otherwise, ¥ k = 1,..., m. 


w i0 " A ij a i,/ p i;3» V k - 1,..*, m. 


For every pair (i,o) e Q, (i,j)) 1 A j.j* 


9 a k £ p m- 


Theorem 3 i n(x) » ji Q + p, x is a lower plane for the NDP, 

4 k 

where \x ^ « 2 r . . .w. * 

k (i,o) e Q 10 

Proof s Consider any pair (i,j) e C. It follows from Proposition 
5 and Proposition 1 that 

A ij x k i^k* x k (irO))# ¥ a k e p ij* (3) 

Summing the inequalities in (3) and losing (l), we get 


E A ii* x k = 2 A n- x *K a ii iPii ^(x, Ci»a)), 

a. eP, , 10 K k=l ^ ^ 

K 1 J 
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which in view of (2) becomes 


m 


2 w ij x k - 7I '( X » o)) • 


k=l 


Multiplying this by r. and summing for all (i,j)eQ, we get 


m 


( 2 w* X, ) < 


r. _ it(X, (i,j)) 


o ± * \ ^ W. Ai j N £-1 ~ -i -I 

(i,j)eQ 10 k=l 10 k (i,j)eQ ia 
Rearranging terms in LHS and using Proposition 2 in RHS yields 


m 


2 ( 
k=l (i,j)eQ 


2 r 13 w k 0 ) x k < L(X - X, Q) - 


or 




o 


+ p 4 x < L (X - X, Q) . 


A procedure, named LP2, to determine h ax) is given 

below. 

The procedure computes the first and second shortest path 
lengths using Mmieka’ s algorithm [2l] . Then, each node pair in 
Q. is considered one by one, p. - is computed and an appropriate 

1 J 

weight is added to each a^ e . 

Proposition 6t The computational complexity of determining 
h 4 (x) is 0(n 3 ). 


Proofs The u. s and v .'s can he determined in 0 (4n 3 ) compu- 
lj iJ 

tations [2l], The shortest path between any pair of nodes 
contains atmost n — 1 arcs, hence 0(n) computations are performed 
in each execution of loop - 1. Since this loop is executed 
n(n-l)/2 times, the complexity of the loop is O(n^). The overall 
complexity of the procedure LP2 is, therefore, 0(n^)* 
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procedure LP2; 
begin 

compute u. ' s and the corresponding g. . r s; 
compute v. ' s; set us* 0; 
for k: = 1 to m do p^s = 0; 
for is = 1 to n-1 do 
for Js = i+1 to n do 
begin 


11 o ! =|i o + r id- u iJ ; 

A ia ! - v id " 'V 

p ij ! ’ °» p « ! = <*> 
i ! = ;)>• 

while q / i do 
begin 

p iJ ! - p io +1 > 

s! - s lq s 


end* 

weights = A la .r 13 /p 1;) ; 
for every e F^ . do p. ^ s = i-i^ + weight 
end {loop ~ 1} 


end; 
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2.4 COMPUTATIONAL RESULTS 

The lower planes derived in this chapter were coded in 
FORTRAN- IV and tested on network design problems of different 
sizes with the aim to assess the sharpness of the lower bounds 
provided and to estimate the computational time requirements. 
Hoang's and Gallo's lower planes were also coded for comparison 
purposes. We also tested the effectiveness of the lower planes ^ 
on a special class of network design problems, known as the 
Optimum Communication Spanning Tree (OCST) [l,13]. In the OCST 
problem, a maximum spanning tree problem is solved to obtain the 
lower bound value after determining p^'s. 

The network design problems were generated randomly with 
assured connectedness. The distances were considered to be 
Euclidean in nature and r. .'s were uniformly generated random 

1 J 

numbers in the range (l, 100). We considered 3 problem sizes 
ranging from 25 nodes and 50 arcs, to 50 nodes and 100 arcs, 
and for each problem size 4 different problems were solved with 
different number of arcs fixed. Out of all the arcs fixed, 
approximately 50 percent were included in the network and the 
rest excluded. 

In Gallo's method as well as in LP1, the lower bound value 
depends upon the order in which arcs are examined for deletion. 
While implementing Gallo's method, we considered the arcs in the 
non- increasing order of the c . values, whereas in IP1 arcs were 

J 
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considered in the nondecreasing order of the number of shortest 
paths between all pairs of nodes using the arc. 

For each problem considered, the lower bounds were 
obtained using Hoang’s, Gallo’s, LP1 and LP2 lower planes. The 
lower bound value is expressed as (lower bound - y. Q ) 100/|i o . 

We also noted the computational times taken by these methods. 
These results are presented in Table 2,1 for the NDP and in 
Table 2.2 for the OCST. 

It is apparent from the tables that LP1 yields lower bound 
values which are generally 1 percent to 3 percent tighter than 
those of Gallo. Even this much improvement is significant in 
the network design algorithms and can curtail the implicit 
enumeration substantially. Both LP1 and Gallo’s method yield 
much sharper bounds than Hoang's method, which is consistent 
with the findings of Gallo [9]. Further, the computational times 
of LP1 are generally 20 percent to 40 percent of those taken by 
Gallo’s method. It may, however, be pointed out that computa- 
tional times taken by Gallo's method can be reduced by incorpo- 
rating reoptimization procedures not used by us. 

The results obtained by LP2 are also encouraging. This 
method takes times which are closer to those taken by Hoang's 
method, but bounds provided are much tighter. Our preliminary 
investigations with the OCST algorithm suggest that use of 
this plane is justified for moderately small sized network 
design problems (upto 25 nodes and 50 arcs). 



Table 2.1s Comparison of various lower planes for the NDP 

(Computational times are in seconds on DEC-1090 system) 
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Table 2.2, Comparison of various lower planes for OCST 

(Computational times are in seconds on DEC 1090 System) 
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CHAPTER. Ill 


EXACT AND HEURISTIC ALGORITHMS FOR OPTIMUM COMMUNICATION 

SPANNING TREE PROBLEM 


3.1 INTRODUCTIONS 

In this chapter, we present the exact and heuristic algo- 
rithms tor the Optimum Communication Spanning Tree (OCST) probler 
Given an undirected network and the amount of communication r^ 
between each pair of nodes i and j, the OCST problem is to deter- 
mine a spanning tree for which the total cost of communication 
between all pairs of nodes is minimum among all spanning trees. 
The cost of communication of a given spanning tree is defined as 
follows s for a pair of nodes i and j, there is a unique path In 
the spanning tree between i and j. The length of the path is the 
sum of lengths of arcs in the path. The cost of communication 
for the pair of nodes i and j is r. . multiplied by the length of 
the path. Summing over all n(n~l)/2 pairs of nodes we have the 
cost of the spanning tree. 

In Section 2, we present the notations used in this chapte 
In Section 3» an exact algorithm based on branch and bound approa< 
is presented. In this algorithm, the O(n^) new lower plane deve- 
loped in Chapter II is used to obtain the lower bounds. In Sec- 
tion 4, we present two-phase heuristic algorithm for solving the 
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OCST problem. In Section 5, another heuristic algorithm, named 
Global heuristic, which is an adaptation of the heuristic algo- 
rithm proposed by Dionne and Florian [5] for the Network Design 
Problem, is given. Dionne and Florian' s algorithm was slightly 
modified for applying it to OCST problem. In Section 6, the 
detailed computational performance of these three algorithms 
along with the networks for testing these algorithms is presented. 


3.2 NOTATIONS'. 


Besides the notations used earlier, the following additional 
notations are used in this chapter. 

S * set of nodes which are included In the spanning sub-graph, 

S * N - S 


R s* 



wtotal 


set of retained pairs 

shortest path distance between the pair i and j, (i e S, 

3 e S), over the spanning sub-graph, 
incumbent objective function value, 
lower bound, 

£_ r. if i e S (amount of _communicat ion of node i with 
jeS .nodes in S) 

E r. . j if i e S (amount of communication of node i with 
jeS nodes in S) 

=* £ £ r. . = total amount of communication over (S, S’) 

ieS' jeS" ^ 


jls v r u lr 1£ 1 £ s 
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I =s set of arcs included in the spanning subgraph, 

E ** set of arcs to be excluded from the spanning subgraph, 

U « A-IUE 

» set of arcs unassigned. 

As already mentioned S is the set of nodes which are 
included in the spanning subgraph and it constitutes a connected 
component with the arcs in the set I. The shortest path dis- 
tances between pairs (i,o), i e S and j £ S over the arcs in 
the sub-graph is given by the u° matrix. 

3.3 THE BRANCH AND BOUND ALGORITHM; 

The branch and bound algorithm is essentially an implicit 
enumeration algorithm. This algorithm proceeds by partitioning 
the set of all solutions into subsets (called sub problems) which 
are of reduced size and are progressively easier to solve. An 
arc is selected and is considered as either included or excluded 
in the spanning tree thereby generating two candidate problems. 
The arc selected is known as the branching arc. A lower bound 
is calculated for the solutions within each subset. The lower 
bound gives the minimum possible value of the objective function 
of a feasible solution in that subset. The best solution avai- 
lable at hand at any stage is called the incumbent solution . If 
for a sub problem, the lower bound exceeds the incumbent value, 

1 that sub problem is eliminated from further consideration. This 
nrocess is known as pruning. If a feasible solution is obtained 
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then the incumbent is updated, if necessary, and further branching 
is stopped for this subset. This is known as fathoming « If the 
sub problem is infeasible then also the subset is eliminated 
from further consideration. The partitioning process continues 
until all subsets are pruned or fathomed. Now, we discuss in 
detail various strategies employed m the branch and bound algo- 
rithm for the OCST problem. Fig. 3.1 contains a flow chart of 
the branch and bound algorithm, a more detailed discussion of 
each of the steps follows. 

Step 1 (Solv ing the S hortest Path Problem) s The all -pair shortest 
path problem is solved using Floyd's algorithm over the graph 
G(N, A). 

Steps 2.3. 13P (Sub p robl em Se lection) s The stack is a list of sub- 
problems yet to be examined by the algorithm. Initially, the 
stack contains a single subproblem which corresponds to the OCST 
problem. New sub problems are added at the top of the stack. 

Finally, when the stack is empty, all sub problems have been 
examined and the incumbent solution is the optimal of OCST problem. 
At the intermediate stages there exist various strategies for 
selecting a new sub problem from the stack. The most appealing 
strategy is the depth first rule whi«h selects the new sub problem 
as the one that was most recently added to the stack. This approach 
results in significant reduction in storage, yields a good fea- 
sible solution early in the search, makes premature termination 
of the algorithm attractive and allows reoptimization of sub problem' 


SOL VE THE 

r pvth Pf :> 8 l ^ 



FIG - 3 .1 FLOW CHART FOR THE BRANCH AND BOUND ALGORITHM 



In view of these advantages we implemented depth first rule in 
the algorithm. 


Steps 4 and 5 (Lower Bounding) ; The lower plane h^(x) developed 
in Section 2.3 is used for computing the lower bound for OCST 
problem. From Section 2.3, we have 

L(X, Q) + p. 3 x < L (X - X, a) 

To get a lower bound, we have to solve the problem 
Min L (X, Q) + x 

such that G(N, X - X) is connected. Equivalently, 


m 


G(N, X - X) is connected, 


Define , 


1 - X, 

> 

Y c X 


if 

It 

o 

l°> 

if 

x k = 1 

l - 

x k 



- 


(3.1) and (3.2) becomes 
m * 

G (N, Y) connected. 
This reduces to. 


m 


Min ( S pf - 


m 


k=l K k=l K K 


G (N, Y) connected. 


(3.1) 

(3.2) 


(3.3) 

(3.4) 
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Thus, we have to find Y c £ which satisfies (3*3) and (3.4). 
This obviously corresponds to a maximum spanning tree (IIST) prob- 
lem, which can bo solved using Kruskal's u 15j algorithm. 

3 

A maximum spanning tree problem is solved with u taken 
as edge weights. Let T be the resulting spanning tree. Then 
the lower bound LB is given by 


LB 


m 

I u 
k~l * 


3 

k 


v 



The procedure LOlffiRBOUND which computes the lower bound 

m 


3 

LB is given. The variable totwgt corresponds to Z n : and 

1- & 

~Z JL 

tweigt to Z !->• f , Heap sort method is employed for forming 

a k £ T 

the maximum spanning tree. 


In computing the termp 0 > max. (u a , u ) is used. This 

-n-vJ 

is because for the pairs (i,j), i e S and 3 e S the shortest 
path between i and 3 will be passing through the arcs included 
in the tree, in the final spanning tree. The shortest path 
distance for these pairs is given b}/- the u° matrix and for the 
other pairs u° matrix contains zero. Since we are using u° 
terms in computing n Q , 8. is initialized to max. [u? , u .} - u . 
to take care of the increased value. 


Step^ 6_( Pruning^ Test) ; After determining the lower bound value, 
an attempt is made to prune the subproblem, that is, to remove 
it from further consideration. The pruning test consists of 
checking LB > Z A and Z' is the incumbent value i.e., the objective 



Procedure LOVffiRBOUND; 
begin 

compute p 0 % - 


n 


n 


2 max. (u?_ , u..).r 


1=1 o=l 


10 10 ' 


10 


set p. s = max. {u? , u. }-u. ; ¥ iUL,.., n, ¥ ji= 1,.,., n; 

-l-J -*~d 

kount s= Jl|| totwgt s= Oj 
for every a k e U do 
begin 

pi= i k ; q *= 0 k J 
lf V = d pi then 

begin 

kount s= kount + 1; 

temp i= d^qS d pq *= -, % 1= d p qi 

NlJss 0} N2s= 0; Rs= 0j 

for is = 1 to n do if + temp then Nls=Nl U {i} 

for o s = 1 to n do if = temp + u o _, then N2t= N2 U {3}; 
for every (i,o) e N1 x N2 do 


qo 


if u = u. + temp + u . then 
2-d 'ad 

begin 

Rs= R U {(i,o)}» 

v ia ! = u ic> T ji := v ia>- u ij := d ia : U J1 != u i3 
end; 

for hs= 1 to n do 

for every (1,0) e R do 

if u + u, < u. then 
ih ho 10 

begin 

u id U ih + V" u io 

end; 

costs = 0; weight s= 0; 
for every (i,o) e R do 
begin 

costs = cost + r. . (u -v. . )j 

if u ic _v ia > then 
begin 

weight s = weight + r^. » 

^id s= u id~ v io J p ji Sa p id 
end | 

u. .s= V. . J u . . * = u. . 

10 10 01 10 

end; 
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\ x ^t - weight; totwgt* = totwgt 
gcost^s^ cost; 


V 1= tamp; := d pq 


end 


+ ■weight; 


end; 

find the maximum spanning tree weight, tweigBt with \i as 
edge weights using Kruskal's algor.-thmj 
LB ;= p- t (totwgt -tweigtO 
end; 
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function value of the best feasible solution enumerated so far. 
If the inequality holds, then the sub problem is pruned, else 
this sub problem is partitioned into two new sub problems. 

Step^ 7„ ^Feasibility Test).; The testing for feasibility can be 
easily incorporated m the LO^ffiRBOUND procedure. If only n-1 
arcs are used m the shortest path calculations between all 
pairs of nodes m the sub problem, then the sub problem is 
clearly feasible to the OCST problem as it forms a spanning 
tree. An arc is used m the shortest path calculation when 
the shortest distance between the pair incident to that arc is 
equal to the length of the arc. The variable 'kount' gives 
the number of arcs used m the sub network. When this ’kount' 
equals n-1, then the solution value is compared with the current 
incumbent value and if found better made the new incumbent. 

Steps, 8^ — 9.8 ^Partitipjringj) : After the pruning test fails, the 
sub problem, say P, is partitioned into two new sub problems. 

To accomplish this, we select a separation arc (a, p) to parti- 
tion the solution set of P into subsets. The global cosb of an 

arc is defined as the increase in the I u _.r. . value 

V'i,o)eQ 10 10 

when an arc is deleted from the network. This is computed in 
the procedure LOWERBOUND and stored in array gcost. The branch- 
ing arc (a, p) is selected such that it has the maximum global 
cost and which connects the already formed subgraph. This 
branching arc criteria makes pruning effective when the arc is 
excluded from the network as it has maximum global cost. The 
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procedure BRANCHINGARC which selects the branching arc (a, f3), 
with a e S and p s S is given. 

Procedure BRANC HING ARC $ 
begin 

mcosts = - 
for every e U do 
if mcost < gcost^ then 
begin 

if i k e S and e S then 
begin 

ca- i^i p**= i k ; mcost := goost^ 
end 

else if i^ e S and 8 S then 
begin 

as= ps= i k ; mcost »= gcost^ 
end 

end 

end; 

The branching arc (a, (3) is first considered as included 
in the spanning tree (x^ = l) and then as excluded from the 
spanning tree (x^^ = 0) . The two sub problems of P correspond- 
ing to x a p = 0 and x^ = 1 are added to sub problem stack in the 
given order. The sub problem with x^o = 1 remains at the top 
of the stack. 
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When an arc (a, p) is included in the spanning tree, all 
the arcs [i, p], 1 e S and i a, if present, should not be 
used in the shortest path calculations in the final spanning 
tree as these arcs form cycles if included in the spanning tree* 
Thus, for the connected sub-network shown in Fig* 3.2 when arc 
[a, p] is selected for inclusion the arcs shown by dotted lines 
are to be deleted from the network. These arcs are deleted 
from the network and the shortest path distances are modified 
using the procedure DELETEARCS given in Section 1.2. The u° 
matrix needs to be updated when an arc (a, p) is added. For 
all pairs (i, p), i e S, u^ and uj^ are set equal to u? a + d^. 
Similarly while back tracking the arcs are added to the network 
and u, u° matrices are modified accordingly. 



A numerical example is solved to illustrate the methodo- 
logy of the branch and bound algorithm. The network contains 
7 nodes and 13 arcs. The distance matrix d and the requirement 
matrix r are given. Tho example network is shown in Fig. 3.3* 
The branching tree enumerated is given in Fig. 3,5. Nodes are 
numbered by the order in which they are enumerated. The lower 
bound obtained at each node is written in a rectangle at that 
node. The left branching indicates inclusion of an arc in the 
spanning tree and the right branching indicates the exclusion 
of that arc. The nodes incident to the branching arc is also 
shown. The Z is initialized to the heuristic solution value of 
91004. The seed node is taken as 4. 
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3,4 TWO-PHASE HEURISTIC ALGORITHM; 

The two-phase heuristic algorithm consists of two phases; 

(1) Tree-building phase, and 

(2) Tree- improvement phase. 

The tree-building algorithm starts with a seed node and 
successively builds the tree by including one arc each time. 

The arc selected was such that it yields a very good feasible 
solution. The method includes that arc m the tree, whose one 
end belongs to the tree already formed (i„e,, belongs to S) 
and the other end belongs to S The amount of communication of 
a node i with all nodes in the other set is given by w^. The 
total communication of nodes in set S with nodes in set £3 is 
given by wtotal. Labels (h) are associated with each node as 
given in Section 3,2 and the arc is selected for which 
h^ + wtotal. d^ + h^ is minimum among all arcs (i, j) such that 
i e S and 3 e S. The procedure SELECTARC which finds the arc 
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to be included in the spanning tree is given* This method ter- 
minates when a spanning tree is formed with n-1 arcs. 

In procedure SELECTARC, all the setting operations are 
of O(n^) complexity and the complexity of loop-1 is also O(n^). 
Thus the complexity of procedure SELECTARC is 0(n 2 ). 


In procedure TREEBUILDING the complexity of determining 

~7 

Uj^’s is 0(n"; using Floyd's algorithm. 0(n) operations are 
performed in loop-1. Since SELECTARC has 0(n 2 ) complexity 
and as loop-2 is executed n-1 times, the complexity of TREE- 
BUILDING is 0(n?). The arrays tail and head contain the end 
nodes of the arcs of the spanning tree and the variable tcost 
the associated cost of the tree. 


Tree 



The spanning tree thus constructed by the tree-building 
algorithm is fed to the tree- imp rc cement algorithm which exa- 
mines each tree arc one-by-one anc tries to interchange it 
with an appropriate non-tree arc so that the cost of communica- 
tion decreases. Whenever a profitable interchange is found, 
it is made and the process is repeated with the changed spanning 
tree. The algorithm terminates when no possible interchange 
yields improvement. Each tree-aro [s,t] is considered one-by- 
one and assuming that the arc is deleted, the two sets of nodes 
S and S are formed. Labels are associated with each node as 
defined in Section 3.2. For each arc in the cutset (S, S), the 
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Procedure SELECTARC; 
begin 

Z r. , 
f oeS 10 
compute w . ♦ = \ _ 

^ Z r i v 
OeS 


if 1 e S 
if l £ S 


, ¥ i!= 1,.,,, n; 


wtotals = Z Z r ; 

isS jeS 10 

c Z h ,vP , if i e S 

f J£ S 0 10 


V is = 1, . . . , nj 


V=j 

L w .u . , if i e S 
jeS 3 13 

mcost « = 

for every [i,ol e S x S do 
begin 

costs r. h. i- wtotal , d . + h ; 

1 10 O’ 

if cost < mcost then 

begin 

p'= i; q. s = Os mcost s- cost 
end 

end { loop-1 ] 

end; 

Procedure TREEBUILDING (seed); 
begin 

determine all pair shortest path lengths u^ ! Sj 
set ss= seed; S: = { s j ; S*s = N- - ^ s } J u gs"“ 
for k: = 1 to n~l do 
begin 

select the arc to be added [p,qj using SELECTARCj 
for every i e S do 


begin 


u iq!“ + d -~; 


pi qi 


lq 


end; { loop-1 } 
tail [k]: = p; head [k]s = q; S;=S U {q}; Si- S - {q} 
end; {loop-2) 


n n 


set tcosts= 

end; 


$• Z Z r ,u ; 
2 i=l o=l 
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cost of communication if that arc is exchanged with the arc 
[s,t] is found and the arc [p, q] with minimum cost is obtained. 
If the arc [p,q] decreases the total cost of communication 
then an exchange is made otherwise another tree-arc is selected 
and the process terminates when all tree-arcs are considered 
for exchange and none yields improvement . The algorithm is 
stated below as procedure TREEIMPROVEMEKT . 

In procedure TREEIMPROVEMENT, loop-1 has 0(n) complexity. 

p 

The setting operations require 0 (n ) computations and loop-2 

p 

also has 0(n ) complexity. Since loop-3 is performed at most 
(n-l) times for each profitable exchange, m total O(n^) compu- 
tations are performed in loop-3. Thus TREEIMPROVEMEKFF has O(n^) 
complexity for every profitable exchange and the variable iter 
gives the number of exchanges performed by the procedure. 

Numerical Example; 

The same example that is considered in Section 3.3 is 
taken to illustrate the two-phase heuristic algorithm. For the 
tree-building algorithm, the starting node is taken as 4. The 
value of w, wtotal, label and the branching arc (p,q) selected 
at each stage is shown in Table 3.1. The tree-improvement 
algorithm performed only one iteration. The tree-arc (4,3) 
is exchanged with the arc (2,3). The solution value obtained 
by the two-phase heuristic algorithm is 91004. In this case 
this is also an optimum solution. The solution obtained by the 
heuristic algorithm is shown in Fig. 3.4. 
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Procedure TREEB1PR0VEMENT j 
begin 

iters « 0 ; ks = Of 

100 i iters = iter+l$ count « - 0 ; 
repeat 

counts = count + 1 ; 

ks= k mod (n-l); ks= k+ 1 ; 

ss= tail [k]; t: = head [k]; S: = 0; S : = 0 ; 

for is = 1 to n do if u? g + d gt = u° t then Ss = S U {i} 

else SsrS U {i}; {loop— 1} 
f r -M» if i e S 

compute w l{ =| r if i e s * 1< = 1 

l deS 

wtotals= S 2 r..; 

ieS oeS 13 

(I w . u? , if i e S 

h£«=\jeS V iJa 1 , »,,, n; 

£ w .u? . if 1 e 3 
v-jeS a ^ 

pcost: = h g + wtotal. d st + h^j mcost: = <»j 
for every [i,j] z (S x S)fl A do 
begin 

cost i = h wtotal . d +h . ; 
if mcost > cost then 
begin 

p;= i; qs = 0 ; mcost % = cost 
end 

end; {loop- 2 } 

until (mcost < pcost) or (count « n-l); {loop- 3 } 
if mcost < pcost then 
begin 

tcosts= tcost - (pcost - mcost); 
tail [k]s= p; head [k]s = q* 

for every ( 1 , 3 ) z (S x S) do u? u° + d + u 1 


goto 100 
end 


10 


ip 


pq 


0 . 

qo' 


end; 
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3.5 GLOBAL HEURISTIC ALGORITHMS 


This heuristic is essentially similar to the branch and 
bound algorithm, except that in computing the lower bound, the 
global cost of deletion of arcs is considered as the weight of 
the arc. In this method, each arc is considered one-by-one, 
and the global cost of deletion of that arc is found and that 
cost is associated with the arc. This will not be a valid lower 
bound, as the effect on a pair may be counted more than once. 

As such there is a possibility of elimination of the optimal 
solution, but Dionne and Florian [5] have observed that the 
probability of over looking the optimal solution is very low 
and an effective heuristic results by considering the above 
weights in calculating the lower bound in branch and bound algo- 
rithm. Since it considers the global cost of deletion of an 
arc as the weight associated with that arc, it is named as 
global heuristic. This algorithm is essentially similar to the 
branch and bound algorithm discussed in Section 3.3* except 
that the following modifications are made in the LOWERBOUND 
procedure. 


1 n n 

i) u is computed as -k £ L u . .r. . 

0 i=l j=l ^ 

ii) ^ is set to Oj V i = 1 ,..., n, ¥ 3 = 1 ,..., n. 

iii) is set equal to the global cost of deleting arc a k 

(i.e., = gcost k ). 
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3*6 COMPUTATIONAL INVESTIGATIONS; 

In this section, we present the computational perfor- 
mance of the branch and bound algorithm, the two-phase heuristic 
algorithm and the global heuristic algorithm. All these algori- 
thms are coded in FORTRAN- IV and implemented on DEC 10 system. 

The networks considered for the testing of these algori- 
thms are generated randomly with assured connectedness. Two 
types of networks are considered for the testing of these 
algorithms. In the first type the distances are considered to 
be Euclidean in nature. The x and y coordinates of the nodes 
are generated randomly in a plane. The plane size is taken as 
(200, 200). Then the required number of arcs are added to the 
network by generating the tail and head nodes of the arcs. The 
lengths of arcs are taken as the Euclidean distance between 
its end nodes. Additional arcs are then added to the network, 
if necessary, to make the degree of each node atleast 2, 

In the second type of networks, the distances are gene- 
rated randomly in the range (0, 200). The arcs are generated 
by the same method as given above. 

For both types of networks , the r M ’ s are generated 
randomly in the range (l, 100). 

The tree building algorithm needs a starting seed node. 
Four nodes are considered as starting nodes, each node from 
each quadrant of the plane having maximum total communication 
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with other nodes. The best solution value among these four 
is taken as the solution provided by the two-phase heuristic 
algorithm. 

For the branch and bound algorithm the starting node is 
taken as that node having maximum total amount of communication 
with other nodes. The lower bound value in the branch and 
bound algorithm depends upon the order in which the arcs are 
examined. The arcs are considered in the nondecreasing order 
of the number of shortest paths between all pairs of nodes 
using the arc. This ordering was found to be very effective 
on the number of nodes examined by the branch and bound algo- 
rithm. For both the branch and bound algorithm and the global 
heuristic algorithm the initial incumbent value is taken as 
the solution value provided by the two-phase heuristic algorithm. 

In Table 3.2, we give the computational performance of 
the branch and bound algorithm and the global heuristic algo- 
rithm for Euclidean networks. The number of nodes examined by 
the algorithm and the maximum depth of the branch and bound 
tree are given. In this table 2 corresponds to the solution 
value obtained by the global heuristic algorithm and 2* the 
optimal solution value obtained by the branch and bound algo- 
rithm. 

In Table 3.3, the results of the two-phase heuristic 
algorithm are given for Euclidean networks. The first 10 prob- 
lems correspond to the 10 problems given in Table 3.2. In this 
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Table, 2 ?~ is the solution value obtained by the first phase 
of the heuristic algorithm (tree building algorithm) and 2^ 
corresponds to the second phase (i.e., tree- improvement algo- 
rithm) . The number of iterations taken by the tree— improvement 
algorithm are also given. 

The results for the branch and bound algorithm, the 
global heuristic algorithm and the two-phase heuristic algorithm 
in the case of random networks are given in Tables 3*4 and 

5.5. 

The branch and bound algorithm when applied to Euclidean 
networks is able to solve problems of size 40 nodes and 69 ares 
in about 40 minutes of CPU time. The times taken by the algo- 
rithm and the nodes examined are found to be exponentially 
growing with the problem size. The sparser problem appear to 
be easier to solve than the denser problems as the lower bound- 
ing is tighter in such cases. The performance of the branch 
and bound algorithm is much more encouraging for random networks 
where problems of size 30 nodes and 54 arcs are solved in less 
than 30 seconds. We believe that this dramatic difference in 
performance results due to large number of near-optimum solu- 
tions which makes Euclidean problems inherently difficult. The 
similar observations have been drawn by researchers in solving 
other combinatorial optimization problem, in particular, the 
traveling salesman problem. 



Table 3* 2: Performance of the Branch and Bound Algorithm and the Global 
Heuristic Algorithm for Euclidean Networks. 
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Table 3.3 s Performance of the Two-Phase Heuristic Algorithm for 
Euclidean Networks. 
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For both the branch and bound and global heuristic 
algorithms, the initial incumbent value is taken as the 
solution obtained by the two— phase heuristic algorithm. No 
improvement is made by these algorithms as the initial incum- 
bent solution is itself the optimal solution for all the prob- 
lems solved. All the time is spent in proving the optimality 
of the solution. This shjws the ineffectiveness of the lower 
bounds available for the network design problems. 

We do not observe any significant difference in per- 
formance between the branch and bound algorithm and the global 
heuristic algorithm. Though the computational times taken by 
the global heuristic algorithm are 1-5 times smaller than the 
branch and bound algorithm, the growth function is again found 
to be exponential in nature. Our experience of global heuris- 
tic algorithm is in contrast with the experience of Dionne and 
Florian [5] for the network design problem whose have reportedly 
obtained ' excellent' results by this method. Our experience 
indicates global heuristic algorithm can not be used to solve 
reasonably large sized Euclidean problems. It is worth noting 
that the lower bounds obtained by the global heuristic algorithm 
are upper bounds on the lower bounds of any branch and bound 
algorithms which considers arc deletion one by one. Henoe 
until we are able to consider block deletion of arcs effectively, 
which will be a major theoretical break through, we can not 
expect the branch and bound algorithm attractive for solving 
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larger sized instances of the network design problems and 
OCST problem in particular. Until then, we shall have to 
take the recourse to efficient heuristic algorithms. 

The two-phase* heuristic algorithm developed by us in 

* *> 

this thesis’ is ..found to be very attractive from the point of 
view of accuracy and computational time. A good solution is 
obtained in the tree-building phase, which is then converted 
into an optimal solution m the tree- improvement phase in all 
the twenty problems considered by us. The computational times 
also appear to be very reasonable. For instance, problems of 
size 100 nodes and 1000 arcs are solved m less than 30 seconds 
of CPU time. The algorithm spends almost equal time in both 
the phases. We attribute this encouraging performance of the 
algorithm to the network topology, the structure of the spanning 
tree. The tree-building and tree- improvement methods exist to 
obtain minimum spanning tree. It is therefore reasonable to 
expect that adaptations of such approaches will yield good 
solutions for the OCST problem. 

We also did some investigations to note the effect of 
variations in the size of the plane in which the points are 
generated randomly and the interval from which r ±d are drawn 
randomly. We did not observe any significant change in the 
performance of any of the algorithm. We also applied the 
algorithm on random planar networks and similar results are 

obtained by us. 
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lu passing, we comment on one possible improvement in 
the branch and bound algorithm. We observed the lower bounding 
routine to be the slowest portion in the algorithm where it 
spends about 90% of the time. The routine currently used 
does not exploit the sparseness of the network. We believe 
that the computational times of the algorithm can be drasti- 
cally decreased if an improved method to determine the effect 
of arc deletion is incorporated. This will allow us to solve 
larger sized problem instances. 
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* : o 5 r s < : 0 m* 

P : 7, ,,, 

l=’(-/r'l<( JJ.EO.u) GO TO 200 
3<rziZ3*>r 

5JIi 

I?(«51S4r.GD.O) GO TO 210 



■9 


?l ( ) 


** 

** 

** 


** 

** 


220 

230 


240 

250 


260 


270 


280 

290 

300 


310 

320 

330 


r sL'VST-f-l 

^JCu i'5 r ) a kRI GH1 

^:(L^3n=s 

rov.* i.r=T-Dr fl )srtw«’tcnr 

z 3 ^ r r m a s 

r f ( < 3 i \f r. 3f . <i«i j qj 

Tip to ie J s tat “to tied 

J^|Sf , ='^oesFti 

r*i^i 

00,220 !Zxi v5 

? = f‘5l ] P /(i GO TO 220 

0 2 3 M .• v 7 ( , / ) 

+ l ? ".PC?, 01 ) CIO TO 220 

3DC r J=P 
roci )=o 
zrmr me 

GO n 460 

: ■ < u£ 

Farr toe ^axlnut* Spanning Tree 

■/;flisnnaPr5K?y N '"‘ rivU Gn ro 340 

;tK5r=L^r/2 

rZ IJiPaxrij, 

v ” T f Q 

00"^4,)”f si ,N 


Jf 

K- 


:(J) 


{£, C ‘h l ?'i? <, * sr 2 G7 TO 27 U . 

IFC^JCO.G \ ’ -,.!■>,* 5:3 rri 27f, 

'ILK JJ = 3lIC<l //U 

W2C Jls^RGCK) 

7sK 

30 rj 260 

3 ;j f 1 Jtj's—p 
- '' : ( j) a, tv ip 

'K^'C-iVO) GO TO 2b0 

Isf!K?rf((3 

<=r * n:cfcj 

{[if; ! ; JSPCU.LT.O) GO TO 2W 
GO "0 78 0 

30 T0 300 

GO 1 3 2 2 1 , 

IFCI,€5.fC3 GO TO 3 30 

< < S P 5 rr 1 5 O ( j ) 4 f - S' E, £ ( K ) 

£T \ l ! , ! j?' ( J - ) ib F,.*(K)) GO TO 31 

* * < * ' *\ i < } ™ i 

" > MiKl )=NK 

; 3 « 0 

^K£K J;=F 
f i i U4( <) = « 

> GJJv, "=rcj'/'« v’i i 

‘ “F» vsr»irr,. , >; l . J u) 

V.-C U-'-UCLCS K 
'■vCK i )=^t4GC«.*cS ' ) 
f.'.SJ =(KSr-l 
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** 

** 


340 




350 

4 


360 

37 0 

380 


390 

400 


‘ ) GH I'U 250 

' . • ro^op-rmsr 

2G) = U . O + EfSUDi'J) *[.3 

^0 v T Ip ’’3 IP 520 

For n tie serp 

Psl 4? l‘{ 

m j 

0 PD = XF? t DD) 

IF? ‘ .GJ.f IfTH 3 ) tD 3 /ft 
I £5 ><• f ,0) t^FPifJ G<) TO 3t>0 
rm.sy.m go ro 363 
s-rrpcrjsp 
?=r 

JD( l , 0)=Jft(T ,PP)+ft?0 

mwr"’” 

go n 35 o 

Z'iH flMO^ 

3ele,t« arcs CP,U) , P in st’fP 

IPAlrtS = 0 

tHCP.*2D,I‘fFTi) 3D I'D 4 30 
D°3 = 3*(P,?) 
d(p,D)=Hfih 

MD,?) = i X’U 
‘it = i 
•12 = 831 

3 3 'UO 1 = 1,8 

£FCC JU,P) *-0P3J,-4S.vja,O)J 
11 = 41 +/ 

5 S T 4 ( 11 ) = f 
33 n 400 

IFCOef + UfU^D.^KP,!)) GO I’D 400 
ss‘rui2) = r 


GO TD 390 


H 0 ) 

V 

■ 1 1 


63 

G 3 

«1 


- 0 ' t ! <• 

i T f . I 
1 r’ ( 7 ; 3 ". 

3 3 in ri* 

r = > : tn) 

DD 41) 12=12,4 
0 = 5-1 ’( i 2 ) 
r ^ c 8 s p \ ? . ( : , .i j „ w « 
1FC( >P)+ HD 

p 4 ss n>u 95 ti 
1 Jj)i< -JUKSJsl 
U -TiSf 1 t J 9 S) =*, 
■)?/<! r 1 f 7 . .. n z t (.’* Q 


ra 

TD 


). 


4 ?f) 

420 


■'nss) «0 TO 410 
, j ) J . JE,0CI,L)) GO TD 410 




20 

ti » 

(u,n 

m £ 

P-JG 

4J 

10 

:j. 

0 } >i 

JO 



\* 

JO 

3s ; 

•; j > 

tPi 





i"' 

’3 

t n 



430 

* •> ; 

’’ » ii 

J£ 





I Ft 

1 ? ' 

tPi.E 

0 . 

J) J 



30 

1 \ \ 

J =1 , 

*13*1 R'l 



l SJ 

»DD 

-HD 





f. = u 

no 

HP 





H j 

, fO 

It 

O 
r— * 

■s* 

*ji) 


440 

HO 

*S l 

, *. ) 

1 

s 3 ( u'^ I ) 

. 4 ! 



J TO 4&0 


)~J 15 1 f=i . I?? I 3, 



Digs 7 


45 a 

457 


** 

** 

** 

** 


** 

** 

** 


470 

4 :lU 


470 


500 
51 u 


520 


SSI) 


Ht/.imsj* jJ 1n 4 *‘» 

J(tj,£)=IS JH 

: rr i n: 

:om us: 

^I?S^? J=SE “, 

rfjdospi 

sr 1 • r,u )=i ; 


17 


1 /' 


‘4t>SP JPM-TQP 


sr*:<mp, , k-<> > 

r S ?}?fU? P J* JS *l 

G £ Jm 0 7 ) = £ 

CFCJNzMiF, 141 n> GJ i'll 120 

,us to- u flosnt it ne*cess^ry 
JFSslFS+i 

y «#» | * 

sTjJJ t = 

07 47 J jsrtl,y 
<5 s 4 K ( 1 , J) ^ fj £ i, n 

ifuwmsta!?) so r j sj y 
in?: If <PH -> * J 

7ii * l ’ a:> ?• 7 i > P A R 

03 >10 'l = J ,^1,41 


' One 


#'*»*•' ^ i, if / v, t ; 

nfa'(n=ru>cn 

go mo js 

7ac<tr ^ r*c 


ca ro b ao 

' > * 4 ,v :u 

< = ’ Cli f J ?, 201 

IF( n . GQ TO 640 

S i m 'IV 

$-} J 

K&iM'm* 10 ro sya 

TOC P, 01=0 
i!>o 4 pi=a 
<s;AR:»nj 4 f ?», oj 

r^c^.i'o. i so to 5?o 
OPi'*^ ' ( r{r :gW) GD W 570 

3 (P, » v ’ 

3 <( >, ^ 1 

'lt=0 

42*4*1 

00 550 r 5 = 1 € M 
kic^jc r,Pi+iVo 

1 it 1,0)1 go ro 54a 

J ( x , .' ) s < < 

j c a # r ) s < < 

‘ * — 1 4 i 

n "i r o- r 

n 'o >vj / 

<* = 3?p + 'KO» J 1 

7?C<<,3£, Ho, £1) GO TO 550 
JC’r 0 =<< 


,41121 


540 



** 
♦ * 


560 

570 

5H0 


610 


620 

630 


640 


5/6 
57 y 


ILiJ [•£'<-}• 2) G3 I'J 
Q7 r r j 

D J 5o J 1 1 s.l , V 3 

r=3::i'ff£i) 

j)0„56) I 2 = ^,0 
0- ; • I?) 

„ <K=.H J ,P) » >*.,»(■ ,) , , 

ho, name 

:ov< l'hjh: 
p = G1 'JiUPl 
7*3 i’’3 53 0 
:d m\ws 

3pJ st’e ^rr (pp/^j 


560 




TO 590 


590 
600 * 


ro 

TO 


630 

530 


5ro>:^ =HpH 

tJPURSsD 

«'t=0 

M2=m 

03 530 jsi,» 

ll = illl fPf> ^° Pa, * MS * nC 1 ^^ 5 } GO Tt 

ssru <i )*r 

73 t3 500 

V2 = V2«i fJ(3v> ' ^ *^E,0( PCM)) GO TO to 1)0 

rr -;j-f 

|FC;t|.£0,D3 G3 

^Mi<i n=i,Mi 

i-"" (rt) 
p 2 . ^ S, _ 1 2 = n 2 f M 

\l[\ R‘U ' P4I RS?? + U( PP ' L) > • NK. Ut I , tO ) 

t *. J ' Cl i ■>'» TR*;)ri 
U 1 )'J( J k 1 m,: );,, 

H < , .')*;( T, (l ) 

JCL, D“M,t ) 
r 9 1 ’ i ; J *' 

V ) * Gu T1 630 

DO 523 JMjsPAIRS 

i = x o lost n 

f-M J)i*f J) 

IS > : . = .fC* , <J+uCK»L) 

1P(a 5 M.Gi.OU ,L>) 
f ( f , .0 = r 5 f i 
J ( ' i » s ) - i ^ J ' 

7 1 ' 1 !• 

^ 5 ' 

5 » s 
I /J: 

MM; n 5 , ;t 
<sf, J v<< ))j 
f, i K ( JO) = J 
33 O'UJ 


GO 


TO 620 


7 * 

. J * 

' ( ) J ) S » * 

1 ^ ** 1 

a37 70J=~I1Oi 
< < ) ) 1 .* 


* * ' J k £ 

<s m:x j.!( 


** 


«• * \ j *i \*'r f j 

!M »5(0*J 
POsj'^GMCO 


3 ?) 
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** 

** 


650 


660 


670 

680 

690 


700 

710 


720 

730 


MO arc tpo, o;i) 

DlPP» )■))*>$') 
m j;,ppi= >P 

ij'U ( pp,00J.L.r.7P0) cjo ro &8o 
V2 = 1*1 

}3 ->5 a i = l,s( 
rsi^ocj.ppjtDpo 
[c( j(x , o5).ut:.t i ju^ 

J( l , P) = I 3U1 

ic joa)=isui 

tfl 3 a H 
SR l K U) = t 
10 n 560 


‘01) so TO 660 


jsoi 3 opo+ J tpo.n 

£FC JC^PJ.LK't 3,01) 
J(P?,I)=IStM f 

JUi PP) = I3J1 
V 2 3 v 2 •! 

s e r m ( v2m 
parvus 

I fc n , so. o) gj ro 
j f( lasr.v) gj ro 
n 670 n=i pi 
i = s £ i um 

00 »70 i? = N2,il 

f.= ■ « (r?) 

. ’ • - • <. r ,pp) atupp. 
IFUSJ1.SS.-IU,0JJ 
T( 1, '01 = 1501 
j(Lv I j=rsji 
CO-H’IMUS 
roP="M-i 
so ’ > >r 

2 "w r r au = 

r f. • •> j i j ■ ( f r • 
i 1 1. a * ) 


> TO 660 


SdO 

68 !) 


0) 

GO 


fQ 670 



Algorithm ' J 

“ ' ■'SD,zsj’Ai,rn v's, 


rtRirUi?., ' 

5-^6 " rr 3 X, 


I~( • V>.$,T.O) 


NO. of fcros = ' , 1 6 / 3 x i 
NG Mamaelr = ',I6/3X 

.. .... „ r o , j , ass.,) Vode =t ' , 1 6 /3 X, 

Zstar = ,1J2,' CPU Titte in iser = ' , I8/3X 

Nodes Examined = ', 16, 'Nodes rationed' =• ',I6/3X 

Vo, of Feasible Solutions = , 1 5/3X, ; 

flo. of It,pi ove-ent r. = ', 16, ' laxl num Septa. = ' , lb) * 
7205 !^ , MP,uiril) ) isI.i«itll) 

Sol it) on . ',5- ‘'P1DC15X, 2014)5 


* j r\ v 

Sf3? 


! r/3x, 


*KJ I£(32,7 Ju) 

Tine Limit Exceeded *5 



** 

** 

** 

** 

** 

** 

** 

** 

** 

*♦ 

** 

** 

** 

** 

** 

#* 

*♦ 

** 

*$ 

** 

#* 

** 

** 

** 

** 

** 

* * 
** 
** 
** 
** 
** 
** 
** 
It* 
** 
** 


** 


** 

** 

** 


10 

20 

3D 

40 

50 

60 


^rojr *m tor irkst - . „ 

Bui LDI i^g Algorithm 

^•V.S.Murty, 




rhhit* Data 


'i 

4 

iss*:> 


* 


X n o d r t ant Variables 


jnS^f r °£ nodes 
Jumoer of arcs 
Seel' value 

2 Sl!SiiS«S 1 ' ,»et»or«s 

tor Random i 4 e t w o r k s 


i£ED' 

r ^ r * < 

H V’ , \J ) 

R C 4 r i) 
K*, 4 ) 

Y 

UIHV I) 

u n 
troru 

uk 3 £ L,|( i ) 

4 'adcd ; 

r\Jbi(M) ; 

R1,R2 ! 

Kl,*Z ; 

run •, 
ji 3 5 3*33 s 

D-aSR(i); 


Starting seed node 

nii?I 25 number : 1000000 
Distance matrix 
gS'-jilrejnent matrix 
|??!£?st path distance matrix 
node 0 £ set S 1 

llves 1 ?L n ^»vf 0f set SB4R 
I 'IF ry if tha? 1 ^ 0 ?^ *? set 
ComnunfsLt *f tha last’ node 

the .other s?t "" 

jj»o*l ji5” n to a “SSd?'' er CUtset ts '5«S) 
v & function value 

‘are q 9 o? S t-h« e c^ ead i node of 

»L? ‘ the spanning tree,* 

arr « a if S f hi le c tail , node of 

lange°o£ t R e spanning tree 1 

III! !? r 5 coordinate 
, for ¥ coordinate 
Minimum degree of node 
toil Is required 
wO »t ti is 4 nodes supplied 
;«*r vor< Generator 


riressR jw ff , 

2 ii 

1 ^ t ; Ri ' x i , x5^<'^ i,/>o, ' l ' C5!) ' 5 ^ 

OA Pft H! ^'^sfAR^l d6dd6o. 1 cocorgna v 

£0 n cJo,2qS,4s* ' 

D0DEG,Rl,R2,Xl,X2,n,r2,D f R,3D£R) 

«* » f t * -1 * <•> if .V *.» « <— — — a. _ -. 

* *«.!- '< 1 / 

:.•>«:• I D£ 

*: UP £(2 2, 40) 


•j! ? 2 A f ' $ f R R t S "J , 0 , x , y , 4 T n ? r r 4 , 

U r • n ^ u ' ^ - tJ !# ' L( - 5 ' 0 j ii A fi f 5 0 ) . r ft It, ( 5 0 ) 


uHUU* ifaii 

33 f 3 i D ■ ’ ^ --- r r - ~<r 

.i . » n » Art., m 


’ 15/3 Xf 


- < a * 

v^''E( r i:f5o) r 5?^?^|^1J|,I lstlc Miorltni*) 

*' JR -i 4 T’{ )A, * 4 o. of noses = *, 15 .* Mo. of arcs s 1 3 , , „ 

«irs.t 22 s?f ed ” Ue = ' ,I5 '’ 5 =t*ort GeS«ator'M 5 «ef 

s, lK'.i‘rf3^ * Phase I is followed by Phase n •/) 

«nisH 
DO 350 IZxl 4 
CADL 1 RTIHSI TIMED 

Initialization 

SEED«ORDEfUrZ) 

X*S£ED 


1 


%I4/) 



** 

** 

** 


** 

** 

** 

** 


** 

** 

** 

** 


0 l<CO = r <Pi«y 

rsjix k,v) +i 

70 T = «m(? + l 

£ p t 1 . 0 . xi 3i r j ^ 

[.i‘KCD = T J 

1 = T 

n vj jo 

so m ‘KC r ) = r 'if nr 

id *j r = i , i 
ii ( £ J = 3 

„ ■)) i=< ,m 

90 JCTf U=0(f, J) 

Fu DVD's 41. jori tfim to obtain U 

01 1 11 c=l , v 
DD Ol 1=1,^191 
m r l? pm 33 i’J 130 

id no t ~ r + 1 , g 
mi.SD.Ki 3.1 n no 
rs c.j) 

T)) (;o Tn 100 

J ( r , J ) = i SOI 

ic i, rms.n 

3 00 : DIF I'M-: 

AM i-t arcs on? oy one 
to corn a scanning tree 

03 HO <D.m=l, MMIVl 
J-'f 
[SOUS’D 

U0 no. moon si ro 120 

ik-ivK/J) 
l( j)=4( JJ+RR 
iSiH = rSO'1tRR 
!=Lii<cj) 

33 n no 

120 aK^)=£SiJ4 

'jnruswrjm+isiM 

Select tD? or ancm no arc (P,Q) 
s’orn tie laoels of. noOes 


1 30 

£ = < 

IFn.aM.I9Fm 

J-K 

30 

140 

X s M = ,) 

IF ( J.SO.TMFm 

30 


rsjs=isj'if it j)*j(i 


T-L J C J3 


150 

160 

33 FD 140 
UBilRC r) = ISIM 

i=onxen 

33 ro no 

:om\M=: 

r=v 

rFCi.’SO.iMFrn 

30 

1 70 

180 

j=v 

ISU1SI) 

IFfl.iO.lMFTtf) 

30 


I5'J^ = rSUM+V(J) 
7 = 019X0) 

*J( f 

190 

33 FD 183 
£»4B£[,iQ) = ISlM 
r=ul 4KCI) 
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2 00 

*f 

tf 

ft 


no 

720 


730 

240 


250 


ft 

ft 

ft 


250 


27 0 


280 

? 9 Q 

BOO 

310 

320 

330 

340 

350 


ft 

ft 

ft 


:uii nr*: 


oage 3 


Seleot to? arc (P,p) 

4 11 = 1 ) >000030 
Ksiirrt 

XFJ MP.UFTY) 33 fa 250 

rF? [. <:o. hptx) on ra 240 

M jci ,j).lt.5ci • 


Kill } TO 230 


r ’ D rn 236 

m r *j ~ £ s in 

0 = r 
7=J 
<< = < 

i=lmc n 
11 n 22 > 
k= > 

7sU 4«( J) 

sn u 210 

3 OKI iJ7 

Mote 4, 0 and 4 ror AL 
>*>’> = KP,0l 

crcueo.riprjf) 00 to 270 
rr=rci, o> 
tff I ) = iH n-RR 

iHi >)=ac r,pj+D?o 

ho, n = jc 1 , 0 ) 
r = Li KU ) 

~3 n 2 6!) 

CD'* r t SEL7E 
r4lL£OiHJf)=P 
fEAPUllHIMO 

r^c<<.£^, r jfJv) yslwcy) 

I F C < < , 4 Z . I 4 P T Y) b UK ( KK ) =1.1 4K ( 0 ) 

LI 'KQ)=X 
k = 0 

C DM f 1 ; M fJC 

ZxQ 

03 273 I = l,M4Tfli 
>0 2 70 J=£ + 1 ,0 
Z=ZtR(I, J)fOCl, J) 

?. ALL' Rrr^E(ri' 1 E 2 ) 
f r^£2 = nvf£2-TTMEl 
4Rri'S( 22,300) S5S0,Z,ITME2 
FORMATC JX, ' Seed' node = MS,' 

1 ' Cine in Msec = ',16) 

:-m r ms 

CALL 1 I MPR0H4, ITER, IJ’IME,Z,HEAD,TAIL,D,R,3) 

4RirSC 22,32 0) Z, I TIME, ITER 

F3RMVTC3X,' 3d j. fa. value = ,110, Ti ne in Msec M6< 

1 40. of Iterations = ,15/) 

4RI r‘S.( 23 , 3 30) S3E0 

FPRUTC 3X,' Seed' node * ',15) 

:34mus 

STOP 

Z43 

oro ir an for TREE - IMPROVEMENT Algorithm 

5'j3R3nn4S I MPR3VC 4 ,IT£R, CTIME, Z, HEAP , TAIL , D , R, J) 
I 4 T 53 ER 4,ITER,rriM£,Z,HFAD(5U),TAIL(50), 


doi.fn. value = 


113, 



oage 4 




** 

** 

** 

** 


10 


?0 

30 


** 

** 

** 


40 


60 

60 


70 

80 


i rijv^nr'i^iD Vd 0 ^ [r(50 ' s 0) , 

i u j(5>J3 ,Sf;r'H5oj , c 50> ' 

} .naoo 

; u. u r r r m 2 ( r t «iK n 

r m = 0 v U1MJ 

< = o 

U’*:4=rr2Ru 

01 100 21 J*J r=l jvn mi 

,)= ISM? <) 

^ f "<='»5§ itE? 5 r g„g°e?iS*Sf2 

3P3= XP,3) 

4£=1 
V2= Hi 

53 3 ) r = l,'| 

GO ro 20 

sirMCMi)=i 

n m 

4-?= 12-1 

SET 4( v2 )sf 

3 3 \ffCV!I(S 

mtjsg 

p It 3 t os oon tras arc C PP , qq) 

?o *o n=3,.u 

jsssmci n 

03.J3 (2=42, 4 
T = S£I\j(i2) 

Ri=p(r, j) 

4 f n = v(l ) f RR 

4( Tlswc J)fRR 

VFOrUz^JULfRR 
3 3 5.9 1 1 = 1 , Ml 
r=s2 rvj( £ 1 ) 
rs.n*o 

i=$f ncf2)' n 

ittmmtw*'" 1 - 3 ' 

n 30 n = \j 2 ,M 
[assn(ll) 

1304=0 

03 73 12=42,4 
7=ssrv(r?) 

r I i8Su.cf) = tsui )i,;JCI,J) 
^jifS^^ISjCPJ^nrAL^OPQ + DAOELfQ) 

03 90 £1=1,41 
£=sspjcn) 

03 90 12=42, M 

j*ss,n?Q 2 ) 

<K = Xl, J) 

GO TO 90 

TSJ4 = bA§EL(a)+WrOTAL*KK+LABELCJ) 

VziWikZi' osn 30 T0 90 

PP = £ 

03 = j. 



90 
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:ur/ ( jjt 
13 )“> * =9339 ('• 


1 I'i ](j|* ) j 3 n iio 

:on l in 

2 Z '• r l 'J tJ <; 

ixci-mie (P«0) *riti (Pp.ooi 
ml m Ut* £ f HK«V0,rAlL and u 

7.s7, - 43 35 P 

ra I uf * o=p-> 

*0=09 

3P9- MPP f ?3) 

3 ,1 U 1 riii.il 

r=s->:n( ro 

33 Ui r 2=12, K 
Q?) 

iJ v * J!‘{5 ’ #P|»J + 9fyHj(py r jj 

:u,u /rn-xrh^) 
i r f *s=t i j-sj-nnsj 
5* =: r f hc m 



oage a 


** 
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*+ 

** 

** 

** 


** 


** 


♦ ♦ 


** 


** 


dir 3rid Random 

•aervor< Pner itors and ng? 

2 i-l l lean Jet *<or< r T ?nerator 

5 HOm J2 'ltJin,'»,rSEtfll, JOODEG,Rl,R2,Xl,A2,n#*2,A. 
Pj r i •; . 2 K R 1 £ R 2 \. X 1 , . A 2 , f 1 l . Y 2 , R & 0 M , 3 R 0 E * C 4 ) , R 4 4 X ( 4 } t 

i 

Jl= J-l 

4X = ( U tX2)/2 
*Y“( ft *■ X 2 ) / P 

snnhi;f?r n SJ ri 211 

t ) s-DRMftU,' 1 > *H-n/2 ') 

sr io 1 

" 1 M <’ L J l ~ 

p\rp sshu jusksd) 
n a a r = t,u 
> 3 i j i = 1 ♦ 1 , 1 

\f t, n = t 'Prv- 
\( J# DspiSTY 
<k- I njpn ,iu) 
rc r, P = x< 

30 «(!,!) = « 

u r r r 3 = 3 
10 a(M)p 
S(V«)=5 
3 ( j, y)=: 

03 3) [= 1 , « 

3£3 «22(I)=J 
R3 H = 0 

33 5 ) J = l,J 
30 RS J« = RSJ'HR( I, J) 

&j nn:.i( nsRsui 

„„ 33 no i = i, v 

70 ixsmmt ,xp> 

IY = lR\'l(Yi,Y?> 

I Ff I. 23.1) G3 TJ 90 
33 Bo J=l, r-i 

80 IFfC rx. 20 .X( J) ) .AMO, ( lY.FQ.YC J))) G3 TO 70 

90 : )<m jus 

xcn =ix 

uu atr 

IFflX.GS 

r f c 1 t . . 

_ 3'JA.K r ) =<< 

100 goj r i jus 

<4ft 33 t?.o r = t ,* 

110 J=IRA J(l, J) 

<=r b\j( i , j ) 

IFUCJ,<). JE.ltfFTY) 30 TO 130 

3SGBSSC J)sOeSrES(J) tl 

3E3B2SOO=OEGReS(K)4l 

rx=xc T)-xc<) 

U = YU)-Y(K) 

xx=ix*r x+i y*ty 

XKal^IXCSORTCXX) ) 

M J,<)=K* 


IX) KK=3 
*/) KK =KK 4 1 


REORDER) 



o*?e 2 


*+ 


1 i) 
HO 


** 


150 


160 


J 7 ft 
IPO 


** 


190 


2 Op 


210 


** 

** 

** 

** 


** 


** 


10 


20 


01 1 l ) 1 = 1 . f 

I-TH '>(? iC !) * Gc; H)iDt:c;) Gn in 150 

ja^ysc jjs'OESRKsljjJ! 


4 s A i i 
K = <(l > — X C I) 
c / = /( r >-f( t> 
xx= t < u <f ! v* ry 
<K = Hr£.(H« r CUH 
\ C f t 1 ) = \ < 

n n in 

: o i * r 4 h 

u Ho r = i 
1 ** K = 0 

n m j = i , j 

HC ,> M )( Gl r 3 lb.) 

5 ?;;lrlK.’<f? rM ’ tJn M m l&J 

t K = J 

: ) ir r o 5 

go td ho 

) j * j j j i , 4 

l P ( 1 n < . G £ . P .1 1' i) fj ( J J ) 20 i’ll 170 

IMHOHsa I) 1/0 

U = J 

CVIi (’lit 
23 iH". H 
R 4 U( [ ) = . 1 U 
ORO-CRf I ) = U 
rjMLi(<o=i> 

2'D 4 I U J*: 

0)2H C = i r 3 
4SX = J 

01 0 0.) r=i n 

,;J Tn 

< = J 

comt ms 
<k=.hosr( ri 

TRO-CRU J = 1R0ER(!C) 

1R0ER(<)=<< 

R MA X (K)=R*fAX(X) 

I 'H X. C I ) = 0 

:omvu£ 

RETJRV 

S4D 

RaHon Meteors Generator 

SIHROim.MS 4G2<V,*,ISEl30,WOPl>BG,Rl,R2 f Al,k2,n,U, h 

nrSGSR Am2,Rl ,H2,RSUM,n f Y?,0RDERC4),R4«(4), 

? 05GREE ( 500 ,romf 5 0), A ISO, 50)^^50,50) 

iMFrysiooioou 

41 = 1-1 

IF(2*4.LE.4*N1) G1 TO 20 
rfRUSf 22,10) 

FOR nr( U, ' M > 4(4-1 )/2 ') 

STOP 

ciimms 

CALL. Si TRAK ISESD) 
oo n i=i,4i 


R, ORDER) 



30 

id 


** 


#* 


bO 

60 


70 


** 


** 


** 

** 

** 

** 


ry 

#R2) 


!\( l, r ) sj\r? 
<K = f ? a *j c m 
R( r , n =<< 

J,n = <x 

U I /«.)-' 
a,? j , n = o 
U n i )'=.') 


03 H ( = 1,4 

r is.* 

*S i ,<: I 

> 1 3) f = I M 

RS I« = R3J‘Ur( L,J) 
r jr<\ ji( r ) = rsum 

•> i r ) r = i , * 

!=r<' mu , n 
<=T um,n 


d-l5trrj?;rf) go to to 



<*=nuui 

ft( f , <)=<K 
A{<, Tj=<< 

BO 

3 0 4 t J Jf I £ 


) J HO 1 = 1 

90 

I “'( 0 2GRFFC 

100 

I=iUM(l f '.| 




GK.NJDIirfG) GU 10 UU 


, J) .'.iS.IMFfx ) GO TO 1 00 
OSGR^Cn^GRL'SCn+l 
OPGR^C T)=UFGPF,2CJ) + l 
■1=1+1 

<K=rR/m u ,n2) 

A ( If T )=<‘\ 

AC J,i)=<K 

go n )o 

no zDvruii*: 

>0 l i d 1 = 1 , 4 
*A •<=) 

)1 1 2 f J J 5 1 N j 

r F C VX.GF.f jrAL(JJ) GU fO 120 

4\*=r3r\U( T) 

<= 7 

120 : jarnuis 

3RU2RU)=< 

no nr a ji(<)=if 

.itsr irm 

2M9 

^uni^tiDi to Generate a Uniform Rartdon 
'J u Ti o ? r s In toe racne IA and TS 

1 M T 2 OF R FiHCTJDM I RAM f I ft, IB) 

<=r« + ( r o-is+imAMU) 

CFC <«G T 0 ) i<L=IB 

FRA'lsK 



, 81588 

81589 

book is to be returned on 


ast stamped 





